Correlated two-photon transport in a one-dimensional waveguide side-coupled to a 

nonlinear cavity 
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We investigate the transport properties of two photons inside a one-dimensional waveguide 
side-coupled to a single-mode nonlinear cavity. The cavity is filled with a nonlinear Kerr medium. 
Based on the Laplace transform method, we present an analytic solution for the quantum states of 
the two transmitted and reflected photons, which are initially prepared in a Lorentzian wave packet. 
The solution reveals how quantum correlation between the two photons emerge after the scattering 
by the nonlinear cavity. In particular, we show that the output wave function of the two photons in 
position space can be localized in relative coordinates, which is a feature that might be interpreted 
as a two-photon bound state in this waveguide-ca\ 

PACS numbers: 



I. INTRODUCTION 

Creating quantum correlations among photons has 
been a subject of major interest for studying the founda- 
tions of quantum theory as well as applications in quan- 
tum information science. As direct interactions between 
photons in free space are extremely weak, generation 
of correlated photons generally requires nonlinear me- 
dia. Electromagnetically induced transparency and pho- 
ton blockade are mechanisms that have been exploited 
to achieve strongly interacting photons [TH3]- Recently, 
studies of two-photon scattering from a two-level system 
inside a one-dimensional (ID) waveguide have also re- 
ported various features of photon correlation 0-0] ■ For 
example, Shen and Fan [6| have discovered the existence 
of two-photon bound states, and Roy 0] Las indicated 
an interesting application of the system as a few-photon 
optical diode. We also note that Shi and Sun |8| have em- 
ployed a formal scattering theory to study multi-photon 
transport in a ID waveguide. 

In this paper, we investigate the correlation properties 
of two photons in a ID waveguide that is side-coupled 
to a nonlinear cavity filled with a Kerr medium (Fig. [1} . 
The nonlinear cavity plays the role of a scatterer. It is 
worth noting that such a Kerr nonlinearity has also been 
employed in coupled cavity array systems for studying 
quantum phase transition |9l-ll4j and nonclassical photon 
statistics [l5l - [l8j . Here we will focus on the transport 
properties of two photons determined by the long time 
solution of the Schrodingcr equation, assuming the initial 
photons are in wave packet forms. We will present an an- 
alytic solution based on the Laplace transform method, 
which has been applied to related photon- atom scatter- 
ing problems [l9j . From the two-photon transmission 
and reflection amplitudes, we show how the two scat- 
tered photons can be correlated in frequency and posi- 
tion variables, with the latter revealing photon bunching 
and anti-bunching effects. Our solution also reveals a 
two-photon resonance condition when the incident pho- 
ton energies match the cavity frequency shifted by the 
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FIG. 1: (Color online). Schematic diagram of the physical 
setup. A ID waveguide is coupled to a cavity filled with a 
Kerr-type nonlinear medium. Photons injected from the left- 
hand side of the waveguide are scattered by the nonlinear 
cavity. As a result, photons are reflected or transmitted in 
the waveguide. 



Kerr interaction. The behavior of transmission and re- 
flection near the resonance will be discussed. 



II. PHYSICAL MODEL 

The physical model under investigation consists of an 
infinitely long ID waveguide and a nonlinear cavity lo- 
cated at the origin (Fig. [T]). We consider a single-mode 
field in the cavity, which couples to right- and left- 
propagating fields of the waveguide via the side cou- 
pling [20|, l2l| so that photons can tunnel between the 
waveguide and the nonlinear cavity. The Hamiltonian 
(with h = 1) of the system is given by 

H = uj c a)a^ a) a) da + / dku: k {f k r k + l k l k ) 



-J I dk 
la 



(1) 



Here a and a* are annihilation and creation operators 
associated with the cavity mode with the resonance fre- 
quency w c . The second term in Eq. (|T|) describes the Kerr 
nonlinear interaction with strength U. The Hamiltonian 
of free fields propagating in the waveguide is described by 
the third term, where l k (P k ) and r k (r£) are, respectively, 
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the annihilation (creation) operators for left- and right- 
propagating waves with wave number k and frequency 
oj k . These operators satisfy the commutation relations 



UU = [f k ,fi} = s(k-k'), [i k A, 



o. 



(2) 



Finally, the last term in the Hamiltonian (fT]) repre- 
sents the coupling between the cavity and the waveguide, 
where J is the tunneling strength. 

For convenience, we introduce even- and odd-parity 
modes operators of the waveguide, 

b k = -j=(r k + h)> c k = -^=(f k - l k ), (3) 

so that Hamiltonian (flj can be rewritten as 

H = H {o) + H {e) (4) 

with 



= I dkuj k c{c k , 

U . 
1 



=uj c a)a + —a^a^aa + 



dkuj k b\b k 



+ g dk{a)b k + b[d). 



(5a) 



(5b) 



Here g = v2J is introduced. We see that the interac- 
tion involves only even modes, and photons in the odd 
modes evolve freely in the waveguide. Therefore we shall 
focus on the calculation of the transport properties of the 
photons in even modes. 

uJcd^a + 



In the rotating frame with respect to Hq 
c J °° dkb k b k , the Hamiltonian can be simplified to 

llf = —d^a i dd + \ dkA k b\b k 
2 j 

POO 

+g / dk(tfb k + bla), (6) 



where A k = u k — uj c is the detuning. In this paper the 
dispersion relation for the modes in the waveguide is as- 
sumed to be linear, i.e., u k — v g k, and we will set the 
speed of light in the waveguide as v g = 1 . 



III. SINGLE-PHOTON TRANSPORT 

As a preparation for finding the solution for two- 
photon scattering, we first consider the single-photon 
problem [22T - [26j . Note that the Kerr nonlinearity has 
zero effect for single photon states. The main purpose 
in this section is to present the single-photon transmis- 
sion and reflection coefficients, which will appear in the 
two-photon solution later in the paper. 

In the single-excitation subspace, an arbitrary state 
can be written as 



!*>(*)>= a(*)|l> £ 



dfc&(t)|o> c |i fc >, 



(7) 



where |l) c |0) stands for the state with one photon in the 
cavity and no photon in the waveguide, and |0) c |l/ £ ) de- 
notes the state with a vacuum cavity field and one photon 
in the fcth (even) mode of the waveguide. The time de- 
pendent variables a(t) and (3 k (t) are the respective prob- 
ability amplitudes. 

By the Schrodinger equation i\<fi(t)) = //j e '|<p(i)), we 
have 



ir(t) = -ig / dk/3 k (t), 



p k (t) = -iA k /3 k (t)-tga(t). 



(8a) 
(8b) 



By performing the Laplace transform defined by f(s) 
Jo°° f{t)^ st dt, Eq. © becomes 



sa(s) 



a(0) 
8pk(s)- /3 fc (0) 



-ig / dk~p k {s), (9a) 
Jo 

-iA k j3 k {s)-iga(s), (9b) 



where a(0) and j3 k (0) are the initial values of the proba- 
bility amplitudes. 

Assuming that initially the cavity is in the vacuum 
state and an incident single photon in the waveguide is 
prepared in a wave packet with a Lorentzian spectrum, 
the initial condition reads 



a(0) = 0, (3 k (0) = 



Gi 



(10) 



where S and e are the detuning and spectral width of 
the photon, and G\ = y/ejiv is a normalization constant. 
The choice of fi k (0) in Eq. (|10[) has the advantage that 
analytic solutions can be obtained conveniently. In addi- 
tion, by noting that e — > corresponds to the monochro- 
matic limit, an incident wave packet of a general form can 
be constructed by coherent superpositions of Lorentzian 
wave packets of various frequencies. 
After some calculations, we obtain 



a(s) 
fa{s) 



1 



lixigGx 



s + x S -i( s + e y 

Gi ( 1 



(11a) 



20 



s + iA k V A fe - S + it s+%5-i(s + e) 



(lib) 



Note that in obtaining Eq. (fTTj) . we have made the ap- 



oo s+iAk 



dA k = 7 /2, 



proximation: J °° ji^dk « f° 
where 7 = 2ng 2 . 

Taking the inverse Laplace transform of Eq. (jll[) . in 
the long time limit, jt/2 — > 00 and et — > 00, we have 



a(t — > 00) 
where 



0. 



t k 



p k (t^oo)=i k Pk(0)e- iAkt , (12 



Afe + i7/2' 



(13) 
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Equation (fl2l) shows that the scattering process results 
in a phase shift 9 k for a single photon with wave vector 
k, where the phase shift is defined by exp(i9 k ) = t k . 

In terms of the left- and right-propagation modes, if 
we assume a photon packet is incident onto the cavity 
from the left, then the initial state can be written as 



V2J0 



dk(3 k (0)fl\(d) 

dk(3 k (0)(b{ + cl) 



(14) 



In the long-time limit, the wave function becomes, 



\<p(t -> 00)} = / dk(3 k (0)e^ Akt (t k rl + rJlm, (15) 



where the transmission and reflection amplitudes are de- 
fined as 



A fc 



A fe +i 7 /2' 



r k 



*7/2 



At + if/2 



(16) 



A similar result has been obtained for the case that a 
single photon is scattered by a two-level system in a ID 
waveguide [22], namely, the transmission amplitude t k is 
zero at the exact resonance. This effect was also reported 
in Ref. [2(| for side coupling with a classical field. 



IV. CORRELATED TWO-PHOTON 
TRANSPORT 

A. Equations of motion and solution 

We now turn to the two-photon scattering problem. 
Since the total excitation number operator of the system 
is a conserved quantity, we can restrict the calculation to 
the two-excitation subspace. An arbitrary state in this 
subspace has the form: 

[>00 

|$(t)) = A{t)\2) c \%)+ / dkB k {t)\l) c \l k ) 
Jo 

poo pp 

+ dp dqC p , q (t)\0) c \l p ,l q ), (17) 
Jo Jo 

where |2) c |0) is the state of two photons in the nonlinear 
cavity and no photon in the waveguide, and |l) c |lfc) is 
the state with one photon in the cavity and one photon 
with wave number k in the waveguide. The last term 
represents the state with no photon in the cavity and 
two photons with wave numbers p and q in the waveguide. 
A(t), B k (t), and C Piq (t) denote the respective probability 
amplitudes. 

By the Schrddinger equation, the probability ampli- 



tudes are governed by: 

A (t) = -ill A (t) - iV2g I dkB k (t) . 



(18a) 

B k {t) = -iA k B k {t) - iV2gA (t) -ig dpC p>k (t), 

Jo 

(18b) 

C P . q (i) = -*(Ap + A q )C p . q {t) - ig(B p (t) + B q (i)). 

(18c) 

We assume that the two injected photons are initially 
prepared in a Lorentzian wave packet. The initial condi- 
tion of the system reads, 



A(0)=0, B k (0) = 0, 
C p , q (0) =G 2 ' ' 



^p — Si + ie A q — 82 + it 
1 1 
A q — Si + ie A p — 82 + ie J 



with the normalization constant 

4e 2 



Go 



V2ti 



1 



{Si-5 2 ) 2 + ^ 



-1/2 



(19a) 



(19b) 



(20) 



Here 5j and ej (j = 1,2) are parameters defining the 
detunings and spectral widths of the two photons. Note 
that C p ^ q has been symmetrized in Eq. (|19b|) because of 
the bosonic character of photons. 

We are interested in the asymptotic solution of C p>q (t) 
in the long time limit. After a lengthy calculation (see 
Appendix fA")) . we obtain for t 7 -1 and e _1 , 



Cp,q{t) — (tpt q Cp^ q (0) + B p q ) e 



-i(A p +A q )t 



(21) 



where t p and t q are defined in Eq. ()13[) . The expression 
of B p q is given by 



B„ 



-2C/G 27 2 



(A, 



q) (A q +q)(Ap + Aq-U + tf) 



1 



(Ap + Aq - 61 - S 2 + lie) 

1 



(A p + A q - Si 



ie + i7r 



A q — 5 2 + ie + i 



(22) 



From Eqs. (f2~Tj) and (|22|) . we notice that the term B p q is 
a non-factorizable function of p and q, implying a correla- 
tion between the two output photons. B p q has a numer- 
ator proportional to the strength of the Kerr nonlinearity 
U in the cavity. In the case U = 0, Eq. (|2"Tj) reduces to 
a simple expression C p>q (oo) = i p t q Cp tq (0) exp[— i(A p + 
A q )t], describing two independent scattered photons. 
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B. Two-photon correlation in frequency variables 

Let us express the results in terms of the left- and 
right-propagating modes. Assuming the two photons are 
injected from the left-hand side of the waveguide, then 
the initial wave function can be written as 



IV>(o)> = 



dpdqC p . q (0) rjf 



t*t 



(23) 



According to Eq. (3) and the solution (f2"Tj) . we obtain the 
long-time wave function, up to an overall phase factor 
exp[-i(A p + A q )t], as 



|^(*->oo)> 



oo POO 





OO pOO 



d P dq{c;]//\ + clJli\m 
d P dq(c;l q flil + c l ;Jlfl)\[ 



(24) 



where 



u p,i 


— tptqCp^q 


(0)4 




(25a) 


r ,n 

u p,q 


— TpTqCpq 


(0)- 




(25b) 


u p,q 


— tpT'qCp^q 


(0). 


- l -B 


(25c) 


°P,9 


— ^qTpCp^q 


(0). 


" 2 B p,q- 


(25d) 



Here 



P,1 



and C: 



p.q 



are, respectively, the two-photon 
transmission and two-photon reflection amplitudes, 
which correspond to the processes in which two photons 
with wave numbers p and q are transmitted into the right- 
propagation mode or reflected into the left-propagation 
mode. In addition, C r v l q {C l p r q ) relates to the process 
where the photon with wave number p (g) is transmit- 
ted into the right-propagation mode and the photon with 
wave number q (p) is reflected into the left-propagation 
mode. 

We point out two interesting situations revealing the 
strong correlation of output photons in the frequency do- 
main. The first situation is achieved by injecting two 
identical photons with Si = 8% = and a narrow spectral 
width e <§; 7. This corresponds to the case when the peak 
frequency of the photons coincides with the resonant cav- 
ity. In this case the two photons are mainly reflected and 
uncorrelated [Fig.^a)], but if they are transmitted, they 
are strongly correlated [Fig. HJb)]. This can be seen by 
the fact that t p = t q = at zero detuning, and hence 
the transmission of both photons is dominated by the 
B p q term. In other words, the two-photon transmission 
near 81 = 82 = is almost entirely due to the nonlinear- 
ity in the cavity. Such a pair of transmitted photons is 
frequency correlated with the two-photon transmission 
probability concentrated along the line A p + A q = 
[Fig. |2£b)]. The uncertainty in the frequencies of indi- 
vidual transmitted photons is of the order of 7, whereas 
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FIG. 2: (Color online), (a) and (b) are plots of |7C" i9 | 2 and 
^Cp^l 2 , respectively, when Si = 62 = 0. (c) and (d) are plots 
of |7Cp^,| 2 and \jCp\ q \ 2 , respectively, when Si — 82 — U/2. 
Other parameters are set as U/y = 10 and e/7 = 0.05. 



the uncertainty in the sum of the frequencies of both pho- 
tons is of the order of e. The smaller e, the narrower is 
the distribution. 

The second situation of interest is two-photon reso- 
nance occurring when the sum of energies of the two in- 
cident photons equals to the energy of a cavity containing 
two photons, i.e., S± + 62 — U. In this case the photons 
can jointly enter the cavity. We show in Figs. [2jc) and 
[2£d) an example with 81 = 82 = U/2 7, where the 
frequency correlation appears more effectively in the re- 
flected amplitude C r p T q , since rj rj (j = p,q). This is 
shown in the narrow distribution in Fig. EJd) . The trans- 
mission part [Fig. EJc)], although they carry most of the 
probabilities, are almost uncorrelated. 



C. Two-photon correlation in position variables 

We now discuss the spatial features of the output pho- 
tons. For simplicity, but without loss of generality, we 
consider the monochromatic limit e — > of incident pho- 
tons. The two-photon transmission amplitude projected 
in position space reads (see Appendix |B|) 



(xi,X2\i>rr) ~ — 16n 2 MG2e 



iE(x c — t) 



0(t-X e )<f>rr(x),(26) 



with 



<t>rr{x) = 



ts 1 ts 2 cos(Sx) — 
7 2 



U 



E-U + 17 



(E + iy) 2 - 4<5 2 



-1*1 



(27) 



Here we have defined x c = (xi+X2)/2 and x = x\ — X2 for 
the center-of-mass and relative coordinates respectively, 
and E = 5i + 8 2 and 8 — (Si — 82) /2. We note that 
Eq. (|26l) is a product of the center-of-mass wave function 
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FIG. 3: (Color online). Spatial features of two-photon trans- 
mission at Si = 52 — 0. (a) |0 rr (a;)| 2 as a function of scaled 
relative coordinate 'yx, for various values of the scaled Kerr 
parameter U/'y. (b) |<^ r r(0)| 2 as a function of the scaled Kerr 
parameter U / 7 . 
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FIG. 4: (Color online). Spatial features of two-photon re- 
flection at 5i = 82 = 0. (a) |</!)ii(a;)| 2 as a function of scaled 
relative coordinate 'yx, for various values of the scaled Kerr 
parameter U/'y. (b) |</>n(0)| 2 as a function of the scaled Kerr 
parameter U/j. 



and the relative wave function <f> rr (x), with exp[iE(x c — 
t)\8[t — x c ) describing the center-of-mass motion of the 
two transmitted photons. The second term of 4> rr {x) is 
a function localized around x = with a width 7 -1 . We 
remark that a similar feature was reported in Ref. @ in a 
photon-atom scattering problem, where the exponential 
decaying function is connected to the existence of photon 
bound states. 

To reveal spatial correlations, we take 8\ = 82 = so 
that the first term of 4> rr (x) can be suppressed. This is 
shown in Fig. [21a) for various values of U. Note that 
|0 rr (a;)| 2 is proportional to the joint probability of pho- 
tons with a separation x, therefore the decaying feature 
corresponds to photon bunching. In particular, the joint 
probability of having both transmitted photons at the 
same position increases with increasing U , but it satu- 
rates when U 7 [Fig. [3jb)] . For the two-photon reflec- 
tion amplitude in position space, we carry out a similar 
calculation and obtain, 

{x U x 2 \i>u) « -16ir 2 NG 2 e- lE ^ +t) 6(t+x c )Mx), (28) 
with 

U 



<t>u{x) = r Sl rg 2 cos(Sx) - 



E-U + ij 



1 



{E + i 7 ) 2 - 4<5 2 



>-\x\ 



(29) 



At 8% = #2 = 8 = 0, the second term causes a dip in 
\4>u(x)\ 2 at x — [Fig. B|a)], which is a signature of 
photon antibunching as the reflected photons repel each 



other. As U increases, the joint probability of having 
both reflected photons at the same position decreases 
[Fig. Sfb)], which is in contrast to the transmitted part. 

Finally we describe the effects of two-photon resonance 
around 81+82 = U discussed in the previous section. For 
simplicity we again consider the case 8\ = 82 here. In 
Fig. [5l we illustrate the dependence of the relative two- 
photon wave function on E = 81 + 82- The effect of two- 
photon resonance is most apparent in Fig.[5fa), where the 
reflected two-photon wave function is strongly localized 
around x = when E = U. Away from the resonance, 
the reflected two-photon wave function exhibits an oscil- 
latory pattern in x [Fig.^c)], which is controlled by the 
two-photon detuning E — U. We also plot the transmit- 
ted two-photon wave function in Figs, [^b) and[5{d), in 
which similar oscillatory patterns are observed. 



V. CONCLUSIONS 

In conclusion, we have presented an analytic solution 
of two-photon scattering inside a one-dimensional waveg- 
uide that is side-coupled to a Kerr-type nonlinear cavity. 
The system provides a scheme to realize correlated two- 
photon transport. The Kerr nonlinearity is found to cor- 
relate photons in frequency variables such that A p + A q 
is a constant, which is a constraint because of the energy 
conservation. In position space, we have shown that the 
Kerr nonlinearity can cause the two photons 'stick' to- 
gether with an average separation distance of the order 
of w g 7" 1 . We may interpret the result as a two-photon 
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bound state, because of the exponential decaying shape 
of relative wave function. However, because of the in- 
terference with single photon processes described by the 
first term in Eq. (25), features of photon correlation may 
only be observed efficiently in certain directions. Finally, 
we note that recent studies of a related topic have consid- 
ered using a single atom as a scatterer [a, However, 
in view of recent progress in achieving a giant Kerr non- 
linearity [H-HI , our work suggests that a nonlinear cavity 
may be an alternative way of regarding the correlated 
two-photon transport problem. 



FIG. 5: (Color online). Dependence of spatial features on 
E — 5i + 82- (a) |0n(:r)| 2 and (b) \4> rr (x)\ 2 . Examples at 
particular values of E are shown in (c) and (d). In these 
figures, we use U = IO7 and 5 = 0. 
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Appendix A: Solution of Eq. (| 18|) by the Laplace transform method 



In this appendix, we give a detailed derivation for the solution of Eq. (| 18|) which governs the transport of two photons 
in a waveguide. We use the Laplace transform method to solve these equations. Under the initial condition (|19p . 
Eq. (fT8|) becomes 



(s + iU)A(s) 



r CO 

iV2g / dkB k (s), 
Jo 



(s + iA k )B k {s) = -iV2gA(s) - ig dpC Ptk (s), 



(Ala) 
(Alb) 

[a + i(A p + A q )]C p , q (s) = C M (0) - ig(B p (s) + B q (s)). (Ale) 

Substituting Eqs. (|Ala|) and (jAlcj) into Eq. (jAlbj) . and making use of the initial condition, we obtain the equation 
for the variable B k (s) as 



[A k -i(s + j/2)]B k (s) = 



2ff 2 



9 



U -is 
2ngG 2 



A„ + Afe — is 



1 



B p (s) dA p 
1 



1 



1 



A k + Si - i (s + e) A fe - 82 



Afe + 82 — i {s + e) Afe — Si + ie 



.(A2) 



where we have made the approximation dp s+i ^ + A k ) ~ 7/^- 
The solution of B k (s) Eq. (IA2[) . by inspection, takes the form: 



B k (s) = 



2ngG 2 



1 



1 



1 



1 



with 



Afe - i(s + 7/2) V Afe + 81 - i(s + e)A k -S 2 + ie A k + S 2 - i(s + e)A k -8i + ie 
F k (s) = -£7 



(l + Ffe( s ))(A3) 



1 



1 



1 



Afe +81- i{s + e) Afe - 8 2 + ie A k + 82 - i(s + e) A k - 81 + ie 



U-is-i-y 8 2 +A k -i(s + e)J 8 2 - i (s + e+ f ) S x + 8 2 - i (s + 2c) 



U-is-i'y Si + Afe - £ (s + e) ) 8 X - i (s + e + 2) 81 + 8 2 - i (s + 2e) 



(A4) 
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Then from Eq. (lAlcp we obtain the following expression 



C p ,q(s) 



G, 



s + i(A p + A q ) 

i'y 



V) 



s+l+iA p \s + e + i(A p + S x ) A p -S 2 +ie s + e + i(A p + 5 2 ) A P -5 X + ie 



s+ 2. + iA q \s + e + i(A q + 5 1 )A q -S 2 + ie s + e + i(A q + S 2 ) A q -S 1 + ie 
2 7 2 1 



s + 7 + iU s + 2e + i (Si + S 2 ) \s + e + % + id! s + e + % + iS 2 J \s + \ + iA p s+^ + iA 



s + 2e + i(S 1 +S 2 )s + e + ^ + i5 1 \s + e + i(Si + A p ) s + \ + iA p s + e + i(Si + A q ) s - 



■ iA, 



s + 2e+ i(Si + S 2 ) s + e + J + iS 2 \s + e + i(5 2 + A p ) s + J + iA p s + e + i(5 2 + A q ) s + \ + iA u 



A p — 5i + ie A q — S 2 + ie A q — Si + ie A p — S 2 + ie 



(A5) 



Until now, we have obtained the expression for C P:q (s). Then we can get the expression for the probability amplitude 
C Ptq (t) by performing the inverse Laplace transform of C p ^ q (s). In particular, since we are interested in the output 
state of the two photons, here we present only the long-time solution of C' p ^ q (t — > oo) as 



C p Jt -> oo) = (t p t q C p .q(0) + Bp^e-^+^t, 



(A6) 



where t p and t q have been defined in Eq. (|13p . and the expression for the correlation term is 



B Pt q — 



-2C/G 27 2 



1 



A p + q) (A, + q) (A p + A q -U + irf) (A p + A q - S 1 - S 2 + 2ie) 

1 1 

(A p + A q - Si+ie + q) (A p + A q -S 2 +ie + q) 



(A7) 



Appendix B: Derivation of two-photon output state in position space 



In this appendix, we derive the wave function of the two-photon output state (|24| in position space. For the 
two-photon transmission process, the corresponding wave function in position space can be written as 



OO /"OO 



(xi,x 2 \ip rr ) = I I dpdqC™ q (xi,x 2 \rlr\\$) 
Jo Jo 

/oo />oo 
/ (t P t q C Ptq (0) + B p .q/A)e-^ +A ^e tA ^e tA ^dApdA q + xi O x 2 . (Bl) 
-oo J — oo 

In Eq. (IB1[) . symmetrization of the two photons has been taken into account by introducing (xi , x 2 \rj,r^\$) = 
A^(e ,ApI1 e ,A « 12 + e pX2 e « Xl ). According to the initial condition given in Eq. (|19b|) . we can get the expression 
for the independent transport part as 

/oo />oo 
/ tptqCp, q (0) e-^+^e^^e^i^dApdAg 
-oo J —oo 

= -%Tt 2 MG 2 t Sl ^t S2 ^ lE+2 ^ x ^ cos(Sx)9(t - x c ), (B2) 

where we introduce the center-of-mass coordinator x c = (xi + x 2 )/2, the relative coordinator x = x\ — x 2 , the total 
momentum E = Si + S 2 , and the relative momentum 5 = (Si — 5 2 )/2. 0(x) is the Heaviside step function and ts 1 -i e 
is defined in Eq. (I16[) . Note that here we have taken the approximation exp[7(xi — t)/2] — > under the assumption 
of 7/2 > e. 

According to Eq. (|A7[> . the Fourier transform of the correlation part B pq can be written as 



/oo ^00 o 
/ _^ e -i{A p +A q )t el A p x lel A q x* d ^ d ^ = A% + ( B3 ) 
-OO J —OO ^ 



with 



1 . 2 [°° [°° 1 1 1 1 



A l = --MUG2I 



2 7-oo 7-oo (A, + if/2) (A„ + i 7 /2) (A p + A g - [/ + i 7 ) (A p + A q -S 1 -S 2 + 2k) 

(B4) 



— - — e^^^dApe^^^dA,,, 

(A p + A q - Si + ie + if /I) p 9 ' 



for I = 1,2. The Fourier transform of the correlation part can be obtained as 

8n 2 MG 2 U f 2 (iB+2e)(3Jc _ t) (iB+2e _ 7) M ( } ( } 

(E — U — 2ie + if) (E + if — i2e) 2 - 4S 2 v ; v 7 

According to Eqs. (|B2[) and (|B5|) . the second term in Eq. (|B1[) can be obtained by making the replacement x c — > x c 
and x — > —x. Then 

(x 1 ,x 2 \^rr) = -lQTT 2 MG 2 e (iE+2t){x "' t) e(t - x c )4> rr {x), (B6) 

with 

Using the same method, we can obtain the wave function for the two-photon reflection state, 

(xi,x 2 \^ u ) = I I dpdqC^ q {x u x 2 ^l\$) w -l&^m 2 e-^ E+ ^ x ^e{t + x c )^u{x), (B8) 

with 



Jo 



<ht(x) = r 5l -ier 52 -ieCos(5x) - — — - U , — jj 2 e {tE+ i' 7)|a:| , (B9) 

E — U — 2ie + « 7 (A + « 7 — i2e)^ — 4ir 

where A/" is defined by (a;i,x 2 |f|f||0) = r\f( e - iA r xl e- iA <> X2 + e- iA * x *e- iA « Xl ). 
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